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Abstract 

The spectrum of light scattered from a Bose-Einstein condensate is studied 
in the limit of particle-number conservation. To this end, a description in terms 
of deformed bosons is invoked and this leads to a deviation from the usual predict 
spectrum's shape as soon as the number of particles decreases. 



PACS number (s): 

03.65.Fd (Algebraic methods), 

03.75.Fi (Bose condensation), 

42.50. Ct (Quantum statistical description of interaction of light and matter) 
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The recent achievements of Bose-Einstein condensate (BEC) with a gas of atoms 
confined by a magnetic trap fl[ has stimulated renewed interest in the question as to 
what signatures Bose-Einstein condensation imprints in the spectrum of light scattered 
from atoms in such a condensate 0. |j . 

As well known, to deal with the dynamics of BEC gas the Bogolubov approxima- 
tion in quantum many-body theory |4j is an efficient approach, in which the creation 
and annihiliation operators for condensated atoms are substituted by c-numbers. One 
shortcoming of this method is that the total atomic particle-number may not be con- 
served after the approximation. Or a symmetry may be broken. To remedy this 
default, Gardiner |J suggested a modified Bogolubov approximation by introducing 
phonon operators which conserve the total atomic particle number N and obey the 
bosonic commutation relation in the case of N — > oo. In this sense, this phonon oper- 
ator approach gives an elegant infinite atomic particle-number approximation theory 
for BEC taking into account the conservation of the total atomic number. 

Along this line, the case of finite number of particle has been recentely investi- 
gated ||, and the algebraic method of treating the effects of finite particle number 
in the atomic BEC has been developed. It results a physical and natural realization 
of the quantum group theory |7j in the BEC systems, whose possibility was already 
suggested in ||], thought in a different manner. 

Here, we shall use the deformed algebra to study the response of a condensate 
with finite number of atoms to the laser light and focus our attention on steady-state 
excitation. 

We consider a system of weakly interacting Bose gas in a trap and a classical 
radiation field interacting with these two- level atoms, where b', b denote the creation 
and annihiliation operators for the atoms in the excited state; a\ a, the creation and 
annihiliation operators for the atoms in the ground state. These operators satisfy the 
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usual bosonic commutation relations. The Hamiltonian of the model reads 



H = hmtfb + h \g(t)tfa + g*{t)b 



(1) 



where g(t) is a time-dependent coupling coefficient for the (classical) laser field coupled 
to those two states with level difference hw. Usually, the time dependence of g(t) is 
given by gexp(-iflt), with Q being the frequency of the laser beam. 

Note that, with the above Hamiltonian, the total atomic particle number N = 
Wb + a) a is conserved. In the thermodynamic limit iV — > oo, the Bogolubov approxi- 
mation 0] is usually applied, in which the ladder operators , a of the ground state 
are replaced by a c-number \fW~c, where N c is the number of the initial condensated 
atoms. As a result, Eq. ([[]) becomes the Hamiltonian of a forced harmonic oscillator. 

Moreover, we have to consider the bath of photon modes, beside the classical driving 
field, so that the total Hamiltonian will be [Bfl 



H = hwb ] b + h^N c g(t)tf + h.c. 

+ £ tt k c{c k + hjN c Y,Z(k) \b ] c k + h.c] , (2) 
k k 

where c k represent radiation modes (of frequency Q k ) which constitute the bath and 
£(fc) is the coupling coefficient pertaining to the internal atomic states. 

Now, by eliminating the heat-bath variables, in view of the Markov approximation, 
in the case of Hamiltonian (Q), it is possible to obtain a quantum stochastic differ- 
ential equation that describes the dynamics of the b mode in the Heisenberg picture 

(3) 



d t b(t) = -iAb(t) - igjN c - Tb(t) + V2Tb in (t) 



where A = vj — Q, and T is the damping rate. Roughly, the latter is given by T = r ]\fW c 
P|, where 7 is the one-atom linewidth M. Finally, bi n (t) is the vacuum noise operator 



(&I n (*)fein(0> = (M*)MO) = o 

(b m (t)b\ n (t')} = 5(t-t'). 



(4) 



The solution of Eq. (^) is well known |J, and in the steady-state regime it becomes 



v2r 

where the semiclassical approximation b(t) — (5 + 5b(t) has been used. In (f|), 5b(u) is 
the Fourier component of the operator Sb(t). 

The spectrum of the light scattered from the atoms is given by the correlation of 
the operators 6'(t) and b(t) ||. Hence, in the steady state, the spectrum of fluctuations 
(5tf (u)5b(u')) results zero everywhere, by virtue of @ and This means that in 
the long time limit, only the equal time correlations survive. 

Let us now come back to the Bogolubov approximation [|J. It destroyes the sym- 
metry of Hamiltonian ([I]), i.e., the conservation of the total particle number is violated 
because [N, if] ^ 0. Then, to preserve the property of the initial model, it is possible 
to determine the following phonon operators || 

B = -$=a!b, B^ = ^=ab ] . (7) 



These operators obey a deformed algebra 0. In fact, a straightforward calculation 
leads to the following commutation relation 

1-2^6, (8) 

where we have introduced a small operator parameter r] — 1/N, which for sufficientely 
large number of atoms is considered as c-number. The algebra defined by Eq. (|D belongs 
to the /-deformed algebra |H|, where in general the deformed operator is related to 
the undeformed one through an operator valued function / as 

B = 6/(6+6) . (9) 

In our particular case, we have 

/(6t6) = y /l- r?(6t6 - 1) , (10) 



and for small deformation we get 



B « b 



1 - ^ f 6+6 - 1 
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With the above in mind, the total Hamiltonian (fj) should be rewritten as 



H = Tiwb ] b + TiVN[g(t)B ] + h.c. 

+ Vkclck + hVNj2t( k ) \ B] °k + h.c. 



(12) 



k k 

Since now iV is a conserved quantity, we can consider it as a c-number (we suppose 
that it coincide with N c , i.e., all the atoms are initially in the condensate). Essentially, 
Eq. ( |i2|) describes the damped dynamics of a deformed oscillator. This is a rather 



cumbersome problem to deal with, as shown in [0]]. Here, we simplify the treatment 
with the following argumentations: in the second term of r.h.s. of Eq. ( Jl2| ) , the nonlinear 
character of B must be taken into account, since it is evidenciate by the radiation-field 
amplitude g; instead, in the last term of r.h.s. of Eq. flP^) , such nonlinear character 
can be neglected due to the weak-coupling assumption with the heat bath. Hence, the 
resulting effective Hamiltonian, in a frame rotating with the laser frequency, is 

HSfg-q 

ICO -f 0' U I 

(13) 



H eff = hAtfb + hVNg (V + b) - hX-^± {b% 2 + 6 f 2 b 



k k 
where we have expressed B in terms of b by means of Eq. ( |TT| ) . 

The nonlinear quantum stochastic differential equation describing the dynamics 
of the 6-mode is now derived from Eq. flT3"| ) 

d t b{t) = t^ll (b 2 (t) + 2b\t)b{t)) - iAb{t) - igVN - r6(t) + V2Tb in (t) . (14) 



It obviously reduces to linear equation (|3]) as soon as r] — > 0. 

The steady state value of the field is given by the solution of the following equation 

\J~Ngrj 







[(3 2 + 2\[3\ 2 ) - iA/3 - igVN - T(3 . 
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(15) 



Of course, the solution of the above equation will be different from that of Eq. (|5|) (we 
refer to the latter as Poo), but they approach each other as soon as N increases, as can 
be seen in Fig. 1. 

The dynamics of the small fluctuations is given by 

d t 5b(t) = A5b(t) + B5b\t) + VZTb^t) , (16) 

where 

A = -iA - T + i^Ngrj ((3 + (3*) , (17) 

B = iVNgrj/3. (18) 

In this case, the solution takes the form 

Sb{u) = -L {[iu - A*} b in (") , (19) 

where 

E(u) = \A\ 2 -\B\ 2 -to 2 -icu(A + A*). (20) 
Finally, the spectrum, by means of Eqs. (|19"D, ( E2~0"|) and @, reads 

S(u) = [ dJ (5 b\ U )8 b(u')) = . (21) 

It is shown in Fig. 2 as a function of iV (besides uj). It trivially vanishes for 77 — ^ 
(i.e., iV — > 00), otherwise it gives a signature of finite number of particles. More 
precisely, we may see that it shows a central peak (typical of the Lorentian shape) 
by decreasing the number of particles. On the other hand, it would be interesting 
to study the transition from this structure to the characteristic Mollow triplet (12 



of a single trapped atom. Unfortunately, our approximations are no longer valid for 
very few atoms, and one should devise a technique to solve completely the problem or 
investigate it numerically. This is a plan for the future study. 

Summarizing, we have seen that the particle-number conservation in BEC requires 
a deformation of the bosonic field, hence the introduction of nonlinearity [1C], which 



may lead (in the limiting case of small number of particles) to observable effects on a 
probe light field. Beyond the oversimplified model used, we retain the measurement 
of the light spectrum in presence of few condensed atoms a promising experimental 
challenge. On the other hand, the use of a BEC with small number of atoms would be 
the subject of next generation experiments |13|| . 

The same aim could be pursued in elementary particle field as well. In fact, the BEC 



may also describe the final state of pions in high-energy-heavy-ion collisions |fL4] , |T5 
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FIGURE CAPTIONS 

Fig. I. A plot of the quantity \\(3\ — \@oo\\ a s a function of N. Values of the 
parameters are: A = and g = 2.5^. Furthermore, arg[/3] = arg^oo] = tt/2 \/N. 

Fig. 2. The spectrum S as a function of u and N. Values of parameters as in Fig. 1. 
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Fig.l S. Mancini and V. I. Man'ko 
Spectrum of light scattered 




Fig. 2 S. Mancini and V. I. Man'ko 
Spectrum of light scattered 



